Two forms of the action variable for the relativistic harmonic 

oscillator 



M.K.Balasubramanya"^ 

Department of Physical and Environmental Sciences 
Texas A&M University-Corpus Christi 

6300 Ocean Drive 
Corpus Christi, Texas, 78412 U.S./^ 
(Dated: February 6, 2008) 

Abstract 

The frequency of a classical periodic system can be obtained using action variables without 
solving the dynamical equations. We demonstrate the construction of two equivalent forms of 
the action variable for a one dimensional relativistic harmonic oscillator and obtain its energy 
dependent frequency. This analysis of oscillation is compared with the traditional solution of the 
problem which requires the use of hyper geometric series. 
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II. INTRODUCTION 

Classical periodic systems can be analyzed elegantly in terms of their action and angle 
variables which constitute a set of canonically conjugate momenta and coordinates. Ac- 
tion variables are proportional to ^pjdqj, where {qj,Pj) are the system's coordinates and 
canonical momenta. For conservative systems they are constants of motion in the manner of 
angular momentum and energy. The frequencies of periodic systems can be found using the 
functional relationship between the action variables and total mechanical energy without 
requiring a complete solution of the dynamical equations. We first summarize here that the- 
ory and situate action variables within its matrix. In Section II we apply this formalism to 
determine the frequency of a non-relativistic simple harmonic oscillator using two equivalent 
contour integral definitions of the action variable. In Section III we extend this formalism 
to a relativistic simple harmonic oscillator, obtain its frequency in two equivalent series rep- 
resentations and demonstrate that it has the correct non-relativistic limit. In Section IV we 
derive the expression for the period of this relativistic oscillator by direct integration and 
compare it with the one obtained using action variables. 

The time evolution of a classical system is governed by its Hamiltonian H which is a 
function of its coordinates Xi, the conjugate momenta pi and the time t. The dynamics of 
such a system is determined by Hamilton's equations of motion 

dH{xi,Pi,t) . dH{xi,Pi,t) 

Xi = ^ , Pi = K ■ (1) 

dpi dxi 

In the case of a particle of mass m moving in one dimension under the influence of a time 
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independent potential energy function V(x) the Hamiltonian is given by if = ^ + ^(^)- 
Such a Hamiltonian is a constant of the motion and is the total energy E of the system. 
Thus, 

|l + yW = £. (2) 
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Canonical transformations are those transformations of one set of coordinate and momen- 
tum {x,p) to another set {X,P) that preserve the form of Hamilton's equations. One such 
transformation is generated by the generating function Wc{x,P) (the suffix C in this and 
similar variables refers to the " classical" rather than the " quantum" nature of the mechanics 
considered), whose arguments are the "old" coordinate, x, and the "new" momentum, P: 

dWc{x,P) ^ dWc{x,P) 

If this transformation transforms the Hamiltonian into a function only of P, then, using 
©, 

dH(P) 

P = = ^ P{t) = P, a constant, 

uX 

dH(P) 

X = = Vo, a constant, X{t) = Vot + Xq. (4) 

Thus X and P evolve very simply in time. The former progresses linearly in time and 
the latter is a constant. Wc{x,P), which generates a canonical transformation in which 
the transformed Hamiltonian is independent of the new coordinate X, is the Hamilton's 
characteristic function. It is related to Hamilton's principle function Sc through Sc{x, P) = 
Wc{x,P) — Et (which generates another canonical transformation, not considered here) 
and, for the case of time independent Hamiltonians, satisfies the Hamilton- Jacobi equation 
obtained by using ([3]) in ([2]): 



1 fdWc{x,P) 
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2,n\ 3. ' 'n^)-^^n (5) 

The use of this method to solve the dynamical problem involves the following steps: (i) 
Define a suitable new constant momentum P, (ii) Integrate Eq. ([5]) to obtain Wc{x, E{P)), 
(iii) Obtain x{X, P) and p{X, P) using Eq. ([3]), and (iv) Express X and P in terms of the 
initial values xq,po and t. 

One particular form of Hamilton- Jacobi theory is particularly suited for the study of 
periodic motion. If an inspection of the Hamiltonian indicates that the motion is periodic, 
then by a particular choice of the new momentum P we can evaluate the period of motion 
without obtaining a complete solution of the dynamical problem. The new canonically 
conjugate coordinate and momentum are chosen to he X = w , P = Jc with 

</c = ^ <lpcix,E)dx, (6) 



where pc{x, E), from ([2]), is ^y2m[E — V{x)] and the integral in phase space is performed 
over one cycle of the periodic motion. Jc is the classical action variable and w the angle 
variable. Since Jc = Jc{E) we can invert it to obtain E = E{Jc)- From Eq. (jlj) the time 
evolution of the new coordinate is w{t) = ut + wq where the constant "velocity" is 

It can be shown that uj is the angular frequency of this periodic motion^. Thus the math- 
ematical problem of finding the frequency of classical periodic motion for a conservative 
system is reduced to that of performing the integral ([6]), solving for E to get E{Jc), and 
evaluating dE/dJc- This is a simple and elegant method for evaluating the frequency of 
a system known to be periodic. Charles-Eugene Delaunay (1816-1872) invented action and 
angle variables in the course of his study of periodicity of lunar motion^. In the early days 
of quantum mechanics, Sommerfeld, in his treatment of the motion of an electron in the 
hydrogen atom, made use of the action variable, and evaluated it using a contour integral in 
the complex coordinate plane^. We follow his example and obtain the classical relativistic 
harmonic oscillator's frequency by evaluating the action variable using two different suitably 
defined contour integrals. There is an equivalent formalism in contemporary quantum me- 
chanics where a quantum version of the action variable, also evaluated by contour integrals, 
is employed to determine the energy eigenvalues of a bound quantum system. This has 
received attention in the last two decadea^i^*^. 

An equivalent definition of Jc is that it is the contour integral 

■^c = Pc{x,E)dx, (8) 

in the complex x plane over a contour C (specified next), where Pc{x, E) is a complex valued 
function of the complex argument x, and is defined as a suitable branch of 

Pc{x,E) = ^2m [E-V{x)]. (9) 

The turning points Xi and X2, both real, are defined by Pc{xi, E) = Pc{x2, E) = 0. These are 
also the branch points of Pc{x, E) in the complex-x plane. We choose a branch cut connecting 
Xi and X2 along the real axis. pc{x,E) is chosen as that branch of the square root which 
is positive along the bottom of the cut. The counterclockwise contour C surrounds the 
two turning points Xi and X2 and the branch cut connecting them. The integral in ([8]) is 
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performed by enlarging C outward to another contour 7, expanding Pc{x-, E) in a Laurent 
series in an annulus that contains 7 and using Cauchy's residue theorem. 

An equivalent construction of the action variable has the —1^ xcip, E)dp form, where 
the clockwise contour C in the complex-p plane encloses the two turning momenta, pi and 
P2, and the branch cut connecting them along the real axis. This second form of the action 
variable is insufficiently explored in the literature, and there is no account of the frequency 
determination of the relativistic oscillator using the action variable. We will demonstrate the 
use of this alternate form of the action variable for the harmonic oscillator in the relativistic 
case, and for comparison and completeness, in the non-relativistic case also. 



III. ACTION VARIABLE FOR NON-RELATIVISTIC HARMONIC OSCILLA- 
TOR 

The Hamiltonian for the non-relativistic harmonic oscillator with spring constant k is 
H = ^ + We will refer to the angular frequency of this oscillator, -y/^, as ujq. We 
first demonstrate the technique of constructing the classical action variable for this simple 
case using the two methods outlined earlier. We will extend it to the case of the relativistic 
harmonic oscillator in a very similar manner. 



A. Action variable in the <f pdx form 



The two turning points of the oscillator, xi and X2 for energy E, are obtained, using 
by setting pc{xi, E) = pc{x2, E) = 0. As functions of energy, they are 



- Xi = X2 

We write the momentum in the form 



2E 



Pc{x,E) 





^ kx"^ 


2m 
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iv mk X 



X / 



(10) 



(11) 



and extend it analytically in the complex-a; plane, with a branch cut connecting the turning 
points Xi and X2 along the real axis. The Laurent series for pc{x, E) in the region |x| > X2 
in powers of ^ is obtained by using the binomial series for the square root: 
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Pc{x,E) = 



^ ^-^ = ivmk X 



2\x) 8 V X 



(12) 



Deforming the contour C outward into a circular contour 7 centered at the origin with radius 
greater than X2, and evaluating the contour integral using Cauchy's residue theorem, we get 

Jc = ^ i^fmk (02 = -\x2^] = —E. (13) 

The angular frequency, from Eq. ([7]), is 01 uj = ujq. For other kinds of oscillators where 
the relation Jq = Jc{E) cannot be inverted to obtain E = E{Jc) in a closed form, we 
can use ^ = The non-dependence of u on the simple harmonic oscillator's energy, 

and thus on its amplitude of oscillation, arises from the purely linear relation between E 
and Jc, which is characteristic of this simple classical periodic system, with its quadratic 
potential energy function. As we will see in the next section there is a richer structure to 
the relationship between Jc and E for a relativistic harmonic oscillator. 



B. Action variable's — <f xdp form 



We demonstrate an alternate, but equivalent, form of the action variable for the simple 
harmonic oscillator. The turning momenta, pi and p2, of the oscillator are defined, using 
(ED, by xciPi, E) = xciP2, E) = 0, or 



-Pl=P2 



V2mE. 



(14) 



We write the coordinate in powers of ^ as 



xcip,E) 



k\ 2m Jmk 



P 



1-1^ 

P 



(15) 



Using (ITSl) we extend xc{Pi E) into the complex-p plane with a branch cut connecting pi 
and P2 along the real axis. We choose that branch of xc{p,E) that is positive just above 
the branch cut. The alternate form of Jc{E) is 



Jc = --^ I xc{p,E)dp 
2vr Jc' 



(16) 



where the clockwise contour C encloses the branch cut. For \p\ > p2, we expand xc{p,E) 
in the Laurent series 



xc{p,E) = Y,a'^p'-'^ 



—I 



'mk 



P 



1 _ i f _ i 1^^^ ^ 

P 



2 \p 



(17) 



Deforming the contour C outward into a circular contour 7' centered at the origin with 
radius greater than p2 and evaluating the integral in f|T6l) . we get 

Jc = i2m) ^ (4 = -]-pA = (18) 

The expected result is that this Jc{E) yields the angular frequency u = dE/dJc = ojq. 



IV. RELATIVISTIC HARMONIC OSCILLATOR 



The relativistic motion of the harmonic oscillator is governed by the Hamiltonian 
H{x,pcr) = \/vc¥?<^ + vn^c^ + \kx^ . The additional suffix R indicates that the canon- 
ical momentum is relativistic. The total mechanical energy of the relativistic oscillator will 
be referred to as and E = E — mc^ is its mechanical energy in excess of its rest mass en- 
ergy. The dimensionless energy related parameter we will use is e = The non-relativistic 
case is characterized by e << 1. 



A. <f pdx form of action variable 



The relativistic orbit equation, obtained from H{x,pcr) = E + mc^, is 



PCR 



\ 



2m 



E - -A;x2 
2 



1 + 



E - Ikx^ 



(19) 



There are four branch points, XjR, j = 1,2,3,4, in the complex-x plane where pcr{x,E) 
vanishes. Two of these are the physical turning points xm and X2r given by 



— XlR — X2R 




(20) 



Their locations in the complex-x plane are similar to those of the turning points Xi and X2 
in the non-relativistic case, given by Eq. (fTOll . The other two branch points of Pcr{x, E) are 
on the real axis at 



X3R = X4R 




(21) 



Their form indicates that these two branch points, unlike xir and X2r, are entirely relativistic 
in character. It is clear that X2r < x^r for all energies. We choose a branch cut of Pcr{x, E) 
connecting xm and X2r along the real axis. Two other branch cuts connect X2,r and 2:4^ 



with X = oo along the real axis. The branch of Pcr{x, E) that we choose for complex values 
of X is positive just below the cut connecting xm and X2r. 
Using Eq. (|T9i) we rewrite pcR in the form 



Pcr{x, E)= ^ a 



3-2j 



mk X 



J=-CO 



X2rY 
X J 



X 



Xar 



(22) 



Comparing this with Eq. ( ITTl) we see that the multiplicative factor (1 + e)5[l — 
which is very nearly 1 for low energies and for |a;| << \x4fi\, modifies the non-relativistic 
Pc{x, E) into the relativistic Pcr{x, E). 
The action variable Jcr{E) is defined as 

1 



CR 



2tt 



PcRix,E)dx, 



(23) 



Cr 



where the counterclockwise contour Cr hugs the branch cut between xir and X2r. Expanding 
PcR in a Laurent series in the annulus X2r < |x| < x^ji we obtain 



Pcr{x, E) = iVmk 



1 + - X 
2 



1 - 



X2R 



^ 1 



X2R 



8 \X4R 



1 



X2R 



64 \X4R 



x^ 



+ powers of x other than x ^. (24) 



The coefficient of in this series, required for evaluating Jcr, is a series in powers of 
(is) ~ 2+^ ^ ^'^^'"^ -'- energies. Deforming the contour Cr outward into 

the circular counterclockwise contour ■jr centered at the origin with its radius less than X4r, 
and evaluating the integral in Eq. (l23l) . we get 

2 



J< 



CR 



E 



1 + 



2 + e 



1 



64 V 2 + e 



(25) 



Comparing Jcr{E) for e << 1 with Jc{E) in Eq. ( fT3l) . we see that the relativistic action 
variable has the correct non-relativistic limit. For low energies, up to order e, we have 
Jcr ~ ;^(1 + |)(1 — jq) ~ (1 + Te^)- Finally, we calculate the relativistic angular 
frequency from the equation 



1 



dJ( 



CR 



1 dJ, 



CR 



dE rnc^ de 



where the energy dependent relativistic factor r]{e) is 

1 



r/(e) 



d 
Te 



1 + 



2 + e 



1 

64 



-r]{e) 



2 + e 



(26) 



(27) 



ri{e) is very nearly 1 for low energies. This expression for uj{i~^, valid for all energies, shows 
a decrease in frequency from the non-relativistic case due to time dilation, for an observer in 
the laboratory reference frame. We also see here the explicit dependence of this frequency on 
the oscillator's energy, and therefore, on its amplitude. This expression is equivalent to the 
series representation of the period of this relativistic oscillator obtained by direct integration 
and shown in Section IV (See Eq. (137|) ). Up to order e, it reduces to ~ ^(1 + |e). 



B. — ^ xdp form of action variable 



Starting with Eq. ([T9|) and extending it into the complex p plane we write the coordinate 



as 



XCR 



2^4 



(28) 



We see that there are two turning momenta pjR, j = 1,2, given by E 
or 

- PlR = P2R = 



PjR'^ + m^c^ = 0, 



2mE I 1 + 2 ) • 



(29) 



These two turning momenta, with a dependence on the relativistic factor ^1 + | invite 
comparison with their non-relativistic counterparts in Eq. ( |T^ . They are also branch 
points of xcr{Pi E) and we choose one of its branch cuts from pm to p2R- The presence of 



\/p^(? + vrfic^ in Eq. (!30l) produces two additional branch points of xcr{p, E) of relativistic 
origin, given by —p^R = p/^R = imc. We choose the second set of branch cuts along the 
imaginary axis connecting each of psR and p^R to p = oo. Further, we choose the branch of 
xcr{p, E) that is positive just above the branch cut connecting piR and p2R. 

The coordinate in Eq. (1281) can be rewritten, using the turning momenta, in a form similar 
to the non-relativitic coordinate Eq. (fT5|) . and suitable for Laurent expansion, as 



xcr{p, E) 



—I 



\mk{l + 



p 



(30) 



where 



f{p\e) 



+ 



Kthr 



1 

21 2 



\P2R 



\P2R 



(31) 
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f{p'^,e) is analytic for |j9^| < |p4i?^|, is 1 at = 0, and thus its Laurent series consists 
only of the non- negative powers of p^. We notice that, in this form, xcr{p, E), apart from 
multiplicative constants, is but a product of p and two varieties of binomial series, one in 
powers of ^ and the other in which converge uniformly in the annulus p2R < p < p^R. 



Expanding the square roots in Eq. (I30l) we get 



XCR = -1 



-2E p 



k P2R 



1 

1 - - 

2 



P2R 
P 



1 f Pm 
P 



1 

16 



Pm 
P 



(32) 

where /j(e), which are inversely proportional to p^jj, are the expansion coefficients in the 
Laurent series of /(p^, e). The coefficient of ^ necessary for evaluating the residue is a power 
series in the parameter ~ (l + |), which, for low energies, is of order e: 



The alternate definition of Jcr{E) is 



271 



JCR = -TT- f XCRiP,E)dp, 



(33) 



where the clockwise contour embraces the branch cut connecting pifi and p2R- We expand 
this contour outward into the circular clockwise contour 7^ centered at the origin with a 
radius less than p4/j. Evaluating the integral in Eq. (|33ll on 7^ using Eq. ( l30l) we get 

E 



Jcr{E) 



1 + 



1 7 

1 eH 

16 256 



128 



(34) 



This expression for the relativistic action variable is a different series representation than 
the one in Eq. fl25l) . and reduces to the non-relativistic Jc in Eq. f|T6|) for e << 1. The 
angular frequency, upt, is given by 

dJcR. 



1 

(^R 



dE 



1 dJcR 
mc^ de 



1 d 

ojq de 



1 



1 + U-T77e 



7 



16 256 



1 

128' 



(35) 



Truncating this series to the first order in e for low energies, we recover the previous result, 



^ ^ -(l + fe). 



V. PERIOD OF RELATIVISTIC HARMONIC OSCILLATOR - TRADITIONAL 
TREATMENT 

We evaluate here the period of the relativistic oscillator by direct integration. The period 
r of a relativistic harmonic oscillator, for all energies, can be obtained in a closed form by 

10 



integrating ^ ^ over one cycle of the motion. For a potential energy function V{x) which 
has the form of a symmetric well and is even in x, 



X2R 



[E - V{x)]dx 



^J[E-V{x)f -m^c'^ 



(36) 



X2R is the relativistic turning point on the right given hy E — mc^ — V{x2r) = 0. For the 
harmonic oscillator, with V{x) = ^kx^, the period integrates to 



r 



mc 



Ve + 2mc2 V2' 2 



37) 



with K — 



E+2mc^ , and the hypergeometric series F^ given by 

ab a{a + l)b{b + l) 



F{{a,b \c\ z) z + 



c(c+l) 



For the weak relativistic case, where E « mc^, we retain terms up to order k^, and 



obtain 



r pa 27r 



e , 1 



2(l + -)5 - 1 + 



1 + -K' 
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Further, |(1 — |) and 



T 



(38) 



(39) 



VI. CONCLUSION 



We have shown the utility of the contour integral definition of the action variable in 
determining the frequency of the relativistic harmonic oscillator. The formalism is easily 
extended to periodic systems in two and three dimensions for separable systems. The non 
relativistic frequency emerges naturally for the e << 1 case. A series representation of 
the frequency of any other relativistic periodic system can be similarly obtained. Further, 
other Hamiltonian models of perodic systems lend themselves to this analysis. The central 
problem in this development is the identification of branch points of the temporally varying 
quantity (e.g., p{x, E) or x{p^ E)) expressed as a function of its conjugate variable and other 
constants of motion. There are four such points each in the two versions of the relativistic 
oscillator considered here, with the expansion parameter for the frequency being the ratio of 
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a "near" and a "far" branch point. The case of a general periodic system in one dimension 
is characterized by its frequency depending on several such ratios of magnitude less than 1, 
with each ratio depending on the system's energy. 
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